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Experiments in coherent nuclear and electron magnetic resonance, and optical spectroscopy corre- 
spond to control of quantum mechanical ensembles, guiding them from initial to final target states 
by unitary transformations. The control inputs (pulse sequences) that accomplish these unitary 
transformations should take as little time as possible so as to minimize the effects of relaxation and 
decoherence and to optimize the sensitivity of the experiments. Here we give efficient syntheses of 
various unitary transformations on Ising spin chains of arbitrary length. The efficient realization of 
the unitary transformations presented here is obtained by computing geodesies on a sphere under a 
special metric. We show that contrary to the conventional belief, it is possible to propagate a spin 
order along an Ising spin chain with coupling strength J (in units of Hz), significantly faster than 
(2J)~^ per step. The methods presented here are expected to be useful for immediate and future 
applications involving control of spin dynamics in coherent spectroscopy and quantum information 
processing. 

PACS numbers: 



I. INTRODUCTION 

According to the postulates of quantum mechanics, the 
evolution of the state of a closed quantum system is uni- 
tary and is governed by the Schroedinger equation. This 
evolution of the state of a quantum system can be con- 
trolled by changing the Hamiltonian of the system. The 
control of quantum systems has important applications 
in physics and chemistry. In particular, the ability to 
steer the state of a quantum system (or an ensemble of 
quantum systems) from a given initial state to a desired 
target state forms the basis of spectroscopic techniques 
such as nuclear magnetic resonance (NMR) and electron 
spin resonance (ESR) spectroscopy [l|, [21 and laser coher- 
ent control [sj and quantum computing [J, . Developing 
a specific set of control laws - pulse sequences- which pro- 
duce a desired unitary evolution of the state has been a 
major thrust in NMR spectroscopy_[l| . For example, in 
the NMR spectroscopy of proteins Q , the transfer of co- 
herence along spin chains is an essential step in a large 
number of key experiments. Spin chain topologies have 
also been proposed as architectures for quantum infor- 
mation processing [t*, 's']. In practice, the transfer time 
should be as short as possible, in order to reduce the loss 
due to relaxation or decoherence. 

The time-optimal synthesis of unitary operators is now 
well understood for coupled two-spin systems [HI, [H, 
[Tsl . [T^ . This problem has also been recently studied in 
the context of Hnear three-spin topologies [H, [13, [11], 
where significant savings in the implementation time of 
trilinear Hamiltonians and synthesis of couplings between 
indirectly coupled qubits were demonstrated over con- 
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ventional methods. In [T^, [13, [3, it was shown that 
the time-optimal synthesis of indirect couplings and of 
trilinear Hamiltonians from linear Ising couplings can be 
reduced to the problem of computing singular geodesies. 
In this paper we extend these methods to linear Ising spin 
chains of size greater than 3. In ^16|, the problem of effi- 
cient propagation of spin coherences was studied. It was 
shown that by encoding the inital state of a spin into cer- 
tain so called soliton operators, it is possible to propagate 
the unknown spin state through a spin chain at a rate of 
(2J)~^ per step. In this paper, we extend these ideas 
to finding efficient ways to synthesize multiple-spin order 
and propagate coherences in linear Ising spin chains. Our 
methods are based on reducing the problem of efficient 
synthesis of the unitary propagator with the available 
Hamiltonians in minimum time to finding geodesies on 
a sphere under a suitable metric detailed subsequently. 
The main result of this paper is that contrary to common 
belief, it is possible to propagate a spin order along an 
Ising spin chain with coupling strength of J (in units of 
Hz) between nearest neighbours, significantly faster than 
(2J)~^ per step. We present exphcit pulse sequences 
for advancing a spin order along an Ising spin at a rate 

~ ^j'^ , which reduces the propagation time to 81% of 
the best known methods. 

Numerous approaches have been proposed and are cur- 
rently used [6| to transfer polarization or coherence 
through chains of coupled spins. Examples are pulse 
trains that create an effective 'XY' Hamiltonian [l^ [2^, 
[2ll | which make it possible to propagate spin waves in 
such chains [23,[H,[23]- In order to achieve the maximum 
possible transfer amplitude, many other approaches have 
been developed, that rely either on a series of selective 
transfer steps between adjacent spins or on simple con- 
catenations of two such selective transfer steps 6, 21]. In 
our recent work, we showed by suitable encodings of the 



2 



spin operators it is possible to transfer the unknown state 
of a spin at the rate of -j time units per step through 
the spin chain . Subsequently, the problem has been 
studied in other contexts, see [23, |2^ and references in 
there. We show that methods presented here are sig- 
nificantly shorter than these state of the art techniques. 
Furthermore, the transfer mechanisms that result from 
the methods presented have a clear geometric interpre- 
tation. 



II. THEORY 

We consider a linear chain of n spins, placed in a static 
external magnetic field in the z direction, with equal Ising 
type couphngs between next neighbors [25l[26| . In a suit- 
ably chosen (multiple) rotating frame, which rotates with 
each spin at its resonant frequency, the free evolution of 
the spin system is given by the coupling Hamiltonian 



SttJ Imzl(m+l)z- 
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FIG. 1: The figure depicts the conventional pulse sequence 
for the creation of multiple spin order as in Equation ([l]). The 
spins are sequentially irradiated with hard {^)y pulses with 
a period oi t = jj between them for the evolution of natural 
coupling. 



We assume that the Larmor frequency of the spins are 
well separated compared to the coupling constant J, so 
that we can selectively rotate each spin at a rate much 
faster than the evolution of the couplings. The goal of 
the pulse designer is to make an appropriate choice of 
the control variables like the amplitude and phase of the 
external RF field to effect a net desired unitary evolution. 

We first consider the problem of synthesizing a uni- 
tary [/, which efficiently transfers a single spin coherence 
represented by the initial operator Ii^ to a multiple spin 
order represented by the operator 2"~^ nm=i ^mz- This 
transfer is routinely performed in multidimensional high 
resolution NMR experiments 0. The techniques devel- 
oped in solving this problem time optimally will be used 
subsequently for other key transfers. 

The conventional strategy for achieving this transfer is 
achieved through the following stages 
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In the first stage of the transfer, operator Iix evolves to 
'2.1iyl2z under the natural coupling Hamiltonian He in -j 
units of time. This operator is then rotated to 2Iiyl2x 
by applying a hard (^)y pulse on the second spin, which 
evolves to '^Iiyl2yl'iz under the natural coupling, followed 
by a hard {■^)y pulse on spin 3 and so on. The final 
state 2"~^(J|^Jj^ Imy)Inz prepared in this manner is lo- 
cally equivalent to YVL^i Imz- The whole transfer 
involves n — 1 evolution steps under the natural Hamil- 
tonian, each taking units of time, resulting in a to- 
tal time of ^^j'-. We now formulate the problem of this 
transfer as a problem of optimal control and derive time 
efficient strategies for achieving this transfer. 



III. THE OPTIMAL CONTROL PROBLEM 

To simplify notation, we introduce the following sym- 
bols for the expectation values of operators that play a 
key part in the transfer. 

xi =< hx > 

X2 =< lllyhz > 
X3 =< 2Ilyl2x > 

X4 =< '^hyhyhz > 



(2) 



X2n-3 —< 2" hyhyhy ' ' ' I(n-l)x > 
X2n-2 —< 2" ^Ilyhylsy ' ' ' I(n-l)ylnz > 



Let X — {xi,X2,X3, . 

system is given by 



, X2r, 



The evolution of the 
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where Uk are the control parameters representing the am- 
plitude of the y pulse on spin fc -|- 1. The problem now 
is to find the optimal Uk(t), steering the system from 
(1, 0, 0, • • • ,0) to (0, 0, • • • , 0, 1) in the minimum time. 

In this picture, the conventional transfer method is de- 
scribed as setting all Uk — for units of time, then 
transfer xi = 1 — s- 2:2 = 1. Now using the control u, we 
rotate X2 = 1 to X3 = 1 in arbitrary small time as control 
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can be performed arbitrarily fast (compared to the evo- 
lution of couplings). Now we set Uk — again, evolve jj 
units of time under natural Hamiltonian, transfer X3 = 1 
to 0:4 = 1 and so on. This way the total time for the 

transfer is ^"^j^"* ■ 

We now show how significantly shorter transfer times 
are achievable, if we relax the constraint that the selec- 
tive rotations on a given spin are only hard pulses. We 
show that if we let the selective operations be carried 
out by soft shaped pulses, along with evolution of the 
coupling Hamiltonian, then we can achieve significantly 
shorter transfer time. The analytical pulse shapes can be 
computed by formulating the problem as computation of 
geodesies on a sphere under a special metric. We divide 
the total transfer into transfer steps. We first transfer 
xi = 1 to Is = a;4 = then transfer X3 = X4 = ^ to 

xc, — xq — ^ and so on and finally at the last transfer 

X2n-3 = X2n-2 = ^ to X2n = 1- The first stcp is de- 
scribed by the following equation. The subsequent steps 
can also be understood by relating them to the following 
equation (we suppress the factor of it J in the front of the 
equation) . 



(4) 
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We make the following change of variables. Let 



ri = Xi,r2 = \/x^ + xj^rs = 
as depicted in the following picture. 

x3 



X4, tanfy — — 

X2 




FIG. 2: The figure depicts the control angle 6 

Using u, we can control the angle 9, so we can think 
of as a control variable. The equation for the evolution 

of ri are 



(5) 



where 



-cosd{t) 
A= ( cos 6l(t) -sm9{t) 
sine{t) 



We now need to find 9{t) for steering the system from 
(1, 0, 0) to (0, in the minimum time. 

We now show that this is equivalent to find the geodesic 
on the sphere from (1, 0, 0) to (0, -^) under the met- 

The time of transfer r can be written as 



/J" y sin^ 9{t) + cos^ 9{t) dt. Substituting for sm9{t) and 
cos9{t) from (O, this reduces to 



dt. 



Thus minimizing r amounts to computing the geodesic 
under the metric 



9 = 



dr\ + drl 



The Euler-Lagrange equations for the geodesic take the 
form ^(^) = 1^ and |r(|t) = Note, rl = 

at ^ ari ' ar\ dt ^ ar^ ' or^ ' ^ 

l — r\ — r1 and along the geodesic, L is constant. We get 



dtrl 

which implies that 



dt' 



(6) 



d rirs - rsri 
dt r\ 



= 



Let / = ''3^1 /I ''3 ^ then / is constant along the optimal 

trajectory. For / to be finite, we need sin6'(0) = 0, i.e., 
0(0) = 0. From Eq. [6l we get 



dt\2^ r2 



Substitute ^ = — coa9{t),^ — s'm9{t), gives 



(7) 



so = /, 



' COS 9{t)) = /sin6'(i). 



9{t) = ft + ip 



Since 9{0) = we get (p — 0. We can now simply search 
numerically for the constant / and it turns out that for 
the transfer from (1,0,0) to (Oj"^i"^); the optimal / 
(in units of ttJ) and the minimum time ti (in units of 
are 
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/ = 0.542, n = 1.94^^4^ 

Now differentiating the expression, f|(i) = ta,n{9{t)) 
gives that u(t) in is constant with a value u = 1.084. 



FIG. 3: The curve in Fig. A, depicts the shortest path con- 
necting the north pole (1,0,0) to a point (0, -^), under 
the metric g. The curve in Fig. B, depicts the shortest path 
connecting (■^, 0) to a point (0, ■^), under the met- 
ric g. 

The subsequent transfer steps, whereby {xk,Xk+i) — 
(72' Tf) transferred to {xk+2,Xk+3) = (72' 72^' 
described by the same equation of the type as @ , except 
we now want to transfer from initial state {-^, '^'^'^^ 
to final state (0, 0, -^) for the intermediate steps and 
finally from (-^, -^,0,0) to (0,0,0,1) for the last step. 
For the intermediate steps, using the same method as 
before, we get / (in the units of nJ) and the minimum 
time T2 (in units of (ttJ)^^) as 

/ = 1.2377, r2 = 1.2692 - (\/2- l)7r. (9) 

The optimal u(t) in Eq. ^ for this transfer is shown in 
Fig. H . 

The optimal transfer for the last step takes the same 
time as in ([8]). The total time (in units of (ttJ)"^ ) is 

-4(n-l)(V2-l)), 

which is approximately 81% of the conventional method 
described earlier. 

We now use these methods for finding efficient ways to 
propagate coherences along Ising spin-chains. 



FIG. 4: Fig. shows the optimal u{t) in Eq. Q for the transfer 
(^,^,0,0) to (0,0,^,^). The time is in the units of 

(^J)-i. 



A. Efficient Transfer of Coherence 

Consider the transfer of coherence 

between the spins 1 and n. 

The fastest conventional method for achieving this 
transfer is based on the concatenated INEPT transfer 
[32,] . whereby bilinear operators '2.1kxl{k+i)z are propa- 
gated through the spin chain using the following steps: 



'2IkxI{k+l)z ^'2IkzI{k+l)x 2/(fc+i)2,/(fe_|_2)z(-10) 

The control Hamiltonian Hi — [Iky + I(_k+i)y) evolves 
for a time corresponding to ^ rotation, (this step takes 
arbitrarily small time), followed by evolution of the cou- 
pling Hamiltonian He for ^ units of time. Under 
this transfer strategy, bilinear operators are propagated 
through the spin chain at a propagation speed of jj. 
We now show how the propagation rate of coherences 
through the chain can be enhanced by preparing suit- 
able spin states and efficient transfer of these states using 
methods described before. Using these ideas, it becomes 
possible to transfer the operators one step along the chain 
in only ^^^^ units of time. 
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FIG. 5: The figure depicts the new pulse sequence for the 
creation of multiple spin order. The spins /s to In-2 are 
sequentially irradiated with shaped pulses as in Fig. 4, with 
pulse duration of T2. The couplings evolve during the pulse 
duration and there is no delay between the sequential pulses. 
The initial rectangular pulse of duration ri on spin 2, creates 
the state {x-i,X4,) = (■^Ji'^j) in (El- The final rectangular 
pulse of duration ri on spin /n-i, creates the desired spin 
order. 



Aife 


'ilkzl{k + l)yl(k + 2)x 


A2fc 


8IkzI(k + l)yjf(k + 2)yjf(k + 3)z 


Asfe 


2I(k + l)xI{k + 2)x 


A4fc 


'il{k + l)xl{k+2)yl(k+3)z 




8IkzI(k + l)yI(k + 2)yI(k+3)x 




I61kzl(k + l)yl(k + 2)yl(k + 3)yl(k+4)z 




^I(k + l)xl(k + 2)yl{k + 'i)x 




^I{k + l)xI{k + 2)yI(k + Z)yI{k+i)z 



TABLE L Relevant operators for transfer of spin order. 



Consider the operator 



A, 



ife 



■A 



2k 



A-: 



3fc 



A, 



ik) 



(11) 



where various operators in the equation are displayed in 
Table I. By application of the Hamiltonian H — He + 
u{t)H2, the operator Afc is advanced one step in the spin 
chain, 



A, 



H 



A 



where H2 = '!rJ{I^k+i)v + I{k+3)y) and u{t) is the 
optimally shaped pulse that transfers (:^):^>OjO) ^ 
(0, 0, in equation ^ in the minimum time. The 



time of this transfer is T2 



and it advances 



the spin operators one step in the spin chain. This is 
best understood by splitting the Hamiltonian H into 
two commuting Hamiltonians Ha — '^Ji'2I(k+2}zI{k+3}z + 



'2I{k+3)zI{k+4}z+uit)I(k+3}y) and Hb = TTj(2IkzI{k+l)z + 

'^I(k+i)zl(k+2)z +u{t)I(k+i)y), whose action can be ana- 
lyzed separately. We identify ((Aife), (A2fc), (-Du), {D2k)) 
and ((Asfc), (A4fc), (A-sfc), {Dik)) with {xi,X2,xs.,Xi). 

Under the action of the Hamiltonian Ha, the evolu- 
tion of the operators {xi,X2TX'i,X4) is exactly identi- 
cal to the one in equation If u{t) is chosen op- 
timally, as described, in Fig. 4, we transfer the state 
((Aifc),(A2fc),(^iO,(^20) - (iiO,0) to (0,0, 
Similarly ((Agfe), (A4fe), (i^sfc), (i?4fe)) = (iiO,0) is 
transferred to (0, 0, 5, ^). So under Ha, A^ is transferred 

to \{D^k + D2k + i'sfc + Dik)- 

Similar analysis shows that under the action of Hf,, 
we transfer ((Du), {D2k), {D^u), {DAk)) = 1, 1, 1) to 
((Ai(fc+i)), (A2(fe+i)), (A3(fe+i)), (A4(fc+i))) = 1,1,1), 
which is the desired operator Afc+i. 

Creation of the operator Ai from the operator Iix can 
also be achieved in many ways. Here we present one 
such preparation, (which is not necessarily time optimal). 
Similarly, the final operator A„_3 can be collapsed on to 
the final state Inx using steps as before. 

(ffc)- (/2„)- (H^)r, 

hx — > '^hyhz — > '^hyhx — > '^hyhyhz 

ihy)^ (ff.)x 1 (12) 

4:Iiyl2yl3x —7=[^Ilyl2yl3x + 8Ilyl2yl3yl4:z)- 

Now by applying a pulse on the first spin, we will 
get a superposition of '^Iizl2yl3x and 8/12 



zhyhyhz- Now, 



letting the natural coupling He evolve for another ^ ro- 
tation period with spin 4 is decoupled, followed by a {T:)y 
pulse on spin 2, we will generate the operator Ai. 

We now show how these ideas can be generalized for 
transferring the unknown state of a spin 1 to spin n. All 
we need to show is that we can simultaneously transfer 
Iix Inx and Iiy Iny, and this ensures that Ii^ is 
transferred to Inz- Furthermore, this ensures that an ar- 
bitrary unknown state of spin 1, {aliz + blix + cliy), is 
transferred to the final state {al^z + binx + ciny)- Con- 
sider the pair of operators (A^, Afc+2), where Afc is defined 
in Eq. pTj) . Then observe, 



(Afc, A 



k+2) 



(Afc.+i, Afc._|_3), 



(k+5)y. 



and 



where H — He -f u{t)TTj{I(^i^^i-jy + If^i^^^-^y 
u{t) is the same optimal pulse shape as in Eq.([9]) and it 
takes units of time to advance the pair (Afc, Afc4.2), 

one step in the spin chain. We can now encode the op- 
erators {Iiy,Iix) via operators (Ai, A3) and transport it 
along the chain at a rate of ^j'^ units of time per step. 
The encoding can be achieved by the following steps: first 
applying a {■^)-x pulse on the first spin, transfer Iiy to 
— /iz, then follow the steps in the last paragraph apply- 
ing the pulses transferring Ii^ to Ai, these pulses will 
also transfer —hz to -^{hy + '^-hxhz)- We then apply 
the pulses transferring 



Ai A2 



A3 
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where Hi = He + u{t)7rJ{l2y + hy),H2 = H^ + 
u{t)T: J {I^y + /sj,) and u{t) is the optimaUy shaped pulse 
that transfers (-^, 0, 0) to (0, 0, -^) in equation 
([5]) in the minimum time. These pulses will transfer 

-^{hy + 2h^l2,) to 



(H2) J_ 



{'^Ilxllyhx + 8Iixl2yl3yl4z) 



(13) 



then by applying a (-1)-;^ pulse on the first spin, we will 
get a superposition of ^Iizl2yhx and Slizhyhyliz and 
let the natural coupling He, evolve for another j rotation 
period with spin 4 and 6 decoupled, followed by a (Tr)^ on 
spin 2, we will generate the operator Ai. Thus, we have 
succeeded in encoding {Iiy,Iix) in (Ai, A3). Compared 
to the recently developed methods, [16| that transport 



the two operators {Iix,Iiy) at a rate of jj per step in 
the spin chain, the time taken by the proposed methods 
is about 81% of these recently developed techniques. 



IV. CONCLUSION 

In this article, we studied the problem of time optimal 
creation of multiple quantum coherences and transfer of 
coherences along linear Ising spin chains. We showed that 
by suitably encoding the initial state of the spin in cer- 
tain operators, it is possible to propagate these operators 
along the spin chain at a rate ~ ; which is signifi- 
cantly faster than the conventional belief that it takes jj 
to advance spin operators in a spin chain. The methods 
developed in this paper may find applications in other 
problems of optimal control of coupled spin dynamics 
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